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^ ■ Abstract 

< ... 

We calculate the three-loop anomalous dimension of the non-singlet transverse operator 

from JV = ltoiV=15. Using some guess we have reconstructed a general form of 
three-loop anomalous dimension for arbitrary Mellin moment N. Obtained result is trans- 
formed into Bjorken-x space by an inverse Mellin transformation. The final expressions 
are presented in both Mellin- N and Bjorken-x space. 
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1 Introduction 



The anomalous dimensions of the non-singlet Wilson twist- 2 operator are known up to 
three-loop order [TJ [2j El HI El El [JJ and the first even moment was calculated at four loops 
[H [9]. There is a similar non-singlet transversity operator, which appear, for example, 
in a study of semi-inclusive deeply-inelastic scattering (SIDIS) [TUl EH E21 E3] and in the 
polarized Drell-Yan process [131 EU E51 ESI EH EH] • The scaling violations of the transversity 
distribution were explored in leading [HI Ell EUJ Ell E21 [23] and next-to- leading order 
[2H [251 [26] . At three-loop order the moments N = 1 to 8 for the anomalous dimension are 
known [271 [281 [291 [30] an d for N = 9 to 14 there are results for 0{nfT F ) contributions [3~T] . 
In spite of the calculations of three-loop anomalous dimension for this operator do not seem 
very complicated with compare to already performed in Refs. [7J E21 E31 EI] it still absent. 
In this paper, however, we want to show how it is possible to obtain such result using 
the methods, which we developed for the computations of the anomalous dimensions for 
composite operators in the maximally extended M = 4 Supersymmetric Yang-Mills (SYM) 
theory [351 ESI EZ1 ESI EH] • Namely, using presented below the first fifteen Mellin moments 
and some guess we have reconstructed a general form of three-loop anomalous dimension 
for the arbitrary Mellin moment N. In spite of the fact that the obtained result seems to be 
correct we should stress, that it is mathematically non-rigorous. More consequent method 
for such calculations was suggested in Ref. [40]. Nevertheless, the result of the recent 
calculations of the next (N = 16) Mellin moment of the three-loop anomalous dimension 
of the non-singlet transversity operator from Ref. [4T] coincides with the prediction, coming 
from our result and can serve as a direct confirmation of the general result presented in 
this paper, at least for even values of N. 

The article is organized as follows: In section [2] we briefly describe the method, which we 
used for calculation of fixed moments of anomalous dimension of the non-singlet transver- 
sity operator. The reconstruction of the three-loop non-singlet transversity anomalous 
dimension from the calculated values is performed in section El The methods of recon- 
struction are described in sections 13.11 and 13.41 The final result is presented in section [4] in 
Mellin- N space. In section [5] we present our result for x-space splitting functions. Finally we 
briefly summarize our findings in section El In Appendix we write down the results for fixed 
moments of anomalous dimension of non-singlet transversity operator from N = 1 to 15. 



2 Calculation 



Our calculations of the fixed moments of anomalous dimension for the non-singlet transver- 
sity operator were performed in the same way, as in [3T] except the last step, where we used 
MINCER 



instead of MATAD 43 



as we interesting only with the divergences. We 
repeat here main futures from [31] and detailed description can be found in original paper. 
The local flavor non-singlet twist- 2 transversity operator is given by 
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(2.1) 



with a^ = (t/2) [Yl v - 7V], A r the Gell-Mann matrices for ST/(3) flavor , V» the covari- 
ant derivative in QCD, q(q) denote the quark and antiquark fields, and the operator S 
symmetrizes the Lorentz indices and subtracts the trace terms. 

Following Ref. [33] we consider the Green's function G^'™' 118 which is obtained by con- 
tracting the matrix element of local operator (12. ip with the source term Jn = A Ml . . . A MiV 

H(p, s)G^ R '™\ r u(p, s) = J N (q t (p) I O^Si,.^ I J (P)> > ( 2 - 2 ) 

where p and s denote the four- vectors of the momentum and spin of the external quark 
line, u(p, s) is the corresponding bi-spinor, A 2 = 0. The contraction with the source term 
Jn = A Ml . . . A MAr allows to write a general expression for the corresponding projector, 
which can be found in Ref. [33]. 

The unrenormalized Green's function has the following Lorentz structure 



fiij TR,ns = ^. (A . p) N-lU a , P S TR,ns (p) + ^ ^ +£ ^ +£ . g ^ ^^M +c 5 A^), (2.3) 



with unphysical constants Ck\k=i...5- To determine E TR ' ns (p) we use, following Ref. [31], the 
following projector (iV c denotes the number of colors) 

£ TR <<» = - \ Nc{Ap)N+1{D _ 2) {n^ P"G^ R n - A.pTr^f^] 

+zA.pTT[a fipP PG i ^ s }} . (2.4) 

Renormalization constants within MS-like schemes do not depend on dimensional pa- 
rameters (masses, momenta) [35] and have the following structure 

Z r (-,a,9 2 )=l + Y,c ( f l \a,g 2 )e-\ (2.5) 



e ' n=l 



where a is the gauge fixing parameter. The renormalization constants define corresponding 
anomalous dimensions: 
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For the calculation of the renormalization constants, following Ref. [3S] (see also Refs. [33 
l4"8l I3H]), we use the multiplicative renormalizability of Green's functions. The renormal- 
ization constants Z? relate the dimensionally regularized one-particle-irreducible Green's 
functions with renormalized one as 

rRenormalized f ~ , «, 9* j = l|m Z T (-,a,g 2 j T BsiIC [Q 2 , Qt B , 0B> 6 ) > ( 2 -7) 

where g B and a B are the bare charge and the bare gauge fixing parameter with 
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re=l 



a B = aZ 3 , Z 3 = Z g . (2i 



To obtain the anomalous dimension of transversity operator we should subtract renormal- 
ization of the external quark field, so 



;pTR,m 



Z„ Z E TR,ns, (2.9) 



where Z q is the renormalization constant for the quark field. 

To find the renormalization constant we compute with the FORM [50] package MIN- 
CER [42J the unrenormalized three- loop E TR ' ns (j>) and determine its renormalization con- 
stant Z S TR, ns from the requirement that the poles in e cancel in the right-hand side of 
Eq. ( 12. 7p . As in our previous calculations of the three-loop renormalization constants 
in supersymmetric Yang-Mills theories [5T] we used a program DIANA [52], which call 
QGRAF [53J to generate all diagrams and the FORM package COLOR [54] for evaluation 
of the color traces. Because the right-hand side of Eq. (12 .1\ contains the bare gauge fixing 
parameter «b, we should perform all calculations up to two loops with the arbitrary gauge 
fixing parameter a (i.e. the propagator of gluon is {g^ v — (1 — oc)qnq v /q 2 )/q 2 ), while for the 
three-loop calculations we used Feynman gauge a = 1. To obtain result for the renormal- 
ization constant from Eq. (12. 7p . we put a = 1 only after expansion in the right-hand side 
of Eq. (12 .7p . The anomalous dimensions 7TRns(-^0 are expanded in powers of a s /(4n) as 



W = £ 7?£.(JV) + (£)%8£.<") + (£)'??£.(*) + ■ ■ ■ ■ (2-10) 



.... 2 
7TR,ns 
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We have calculated the fifteen first moments of the three-loop anomalous dimension 7^ 
of the non-singlet transversity operator (12. ip . which are listed in Appendix. 
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3 Reconstruction 

Having in hand fifteen first values for the three-loop anomalous dimension of the flavour 
non-singlet transversity operator we have tried to reconstruct a general form of anomalous 
dimension equipped with our experience in such exercise, which we have done during 
the computations of higher orders anomalous dimensions of composite operators in the 
maximally extended Af = A SYM theory [351 ESI E3 ESI EH] • 

As well known from the direct calculations in QCD the anomalous dimensions of twist-2 
operators are expressed in terms of harmonic sums |55l [561 EZ], following Ref. [56] 
recursively defined by 

S± m {N) = Y^ - m > S , ±mi,m 2 ,...,m fe (A r ) = ^ - mi -S'ma,... > m fc («) • (3-1) 

4=1 % 1=1 % 



The sum of the absolute values of the indices rrii defines the weight of the harmonic sum. 
In the n-loop anomalous dimensions one can encounter sums up to weight 2n — 1. Thus, 
the harmonic sums form the basis and if we calculate a sufficient number of fixed moments 
of the anomalous dimension we can restore a general form of anomalous dimension. We 



used such method earlier for the computations of anomalous dimensions at two-loop 
three-loop [51], four- loop [37] and five- loop [5S] orders for the twist-2 operators and even 
at six-loop [39J order for the twist-3 operator in M = 4 SYM theory. 

To find a minimal basis we used the maximal transcendentality principle, which was 
discovered during the calculations of the eigenvalue of the kernel of Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) equation [SHI ESI ED] at the next-to-leading logarithm approximation in 
M = 4 SYM theory [61J. It was found, that the final result contains only the most com- 
plicated functions from the full result in QCD [52] ■ Assuming a deep relation between 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [631 El El EH] and BFKL equations, 
the authors of Ref. [61] suggested that the eigenvalues of the anomalous dimension matrix 
of twist-2 operators in J\f = 4 SYM theory should contain the most complicated harmonic 
sums, i.e. the harmonic sums only with the maximal weight (maximal transcendentality). 
This hypothesis was confirmed by direct calculation at the two loops in Ref. [35J. So, 
the result in the maximally extended JV = 4 SYM theory can be obtained from the cor- 
responding result in QCD, if we take from the QCD result only those terms, which have 
the harmonic sums with maximal transcendentality (with maximal weight). Note, that 
the eigenvalues of the anomalous dimension matrix of twist-2 operators in supersymmetric 
theories can be expressed through one function with shifted argument. In Ref. [61J such 
function was called an universal anomalous dimension. The universal anomalous dimen- 
sion of the Af = 4 SYM theory is contained in any diagonal element of the anomalous 
dimension matrix of twist-2 operators in any four-dimensional gauge theory, for example, 
in QCD, but with corresponding colour factors. 

The anomalous dimension of the non-singlet transversity operator should contain the 
same universal anomalous dimension as a part of full result. Moreover, one can see, that 
one- and two-loop anomalous dimensions of the non-singlet transversity operator [191 EDI 
EH HH E21 [2H [251 EE] differ only by few terms with compare to the universal anomalous 
dimension [HH ES]- We assume, that this property will be hold also at three loops. 

Let's see how it works at two loops. We take the result for the two-loop anomalous 
dimension of the non-singlet operator from Eq. (3.5) of Ref. [7] 
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N. 



S 2 + 2S 3 
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Si + 45i 



2 Si a + 25*2,1 + S3 



-4 C A C F ( 2 N + S 3 - H - 25_ 3 - f Si + (N_ + N + ) 
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(3.2) 



7^-(A0= 7 y s j+ (A0 



(1)H 
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) ((N_ - N + ) 



So — S3 



2(N_ + N + -2)S] 



(3.3) 



where N± S^ = Sj^(N ± 1) and make substitutions N_ — > 1 and N + — y 1. As result, we 



obtain 

1 10 „ 2 \ „ o/.„ „ 3 



l£> + (N) = 4C> / (- - —Si + -5 2 ) +AC/(AS. 3 + 35 2 - - + 25 3 - 85i,_ 2 

- 4S 1)2 - 452,1 - 25 3 ) + 4 C A C F (2 5 3 - ^ - 25_ 3 + y 5i + 45 r i,_ 2 - y 5 2 ) , (3.4) 
7^(A0= 7l ( i )+ (A), (3.5) 

which should be compared with the two-loop anomalous dimension of transversity operator, 
calculated in Refs. [2H E3 [26]: 

tSUW = 4 C>V (1 - y 5 X + 2^ + 4 ^2 ^_ 3 + 3S2 _ 3 + ^ _ g ^_ a 
- 45-1,2 - 452,1 - 25 3 ) + 4 C A C F (2 5 3 - ^ - 25_ 3 + y 5i + 4Si,_ 2 " y S ? 

The difference is only in the last term. So, we can expect, that such simplification will 
work at the three-loop order too. 

3.1 C F nf contribution 

We start the reconstruction of the three-loop anomalous dimension of non-singlet transver- 
sity operator from the most simple contribution proportional to C F n?, which already known 
from previous calculations [281166] . Again, we take Eq. (3.7) from arXiv version of Ref. [7] 



7 ^+ JN) = lQC F nf(— - — 5i + -S 2 + (N_ + N 



2 9 11 9 + 1 9 
9 Sl ~ 5T 2 + IT 3 



(3.7) 



In this equation we make substitutions N_ — > 1 and N + — > 1 and obtain 

,(2) ^_1fi^^ 17 l - 5 - l 



7 ^ JN) = 16C F n f 2 [ Si S 2 + -S 3 ). (3. 

'ns,c F n/v 1 F ' \144 27 27 9 / v 



Then we subtract the last equation at iV = 1, . . . , 15 from calculated values Eqs. (lA.lj) - 
(IA.15I) with the following result 

2 r4 4 2 2 4 4114 4 2 2 4 4 1 l 
prif 19 '27 '27 '45' 135' 189' 6^ ' ^ ' 4^ ' i% ' 2^ ' 3^1 ' 819 ' 945 ' 270 J ' ^ ' ' 

To reconstruct a full N dependence for obtained numbers we should write a suitable 
basis from the harmonic sums and solve the system of equations to find coefficients in the 
front of sums in ansatz. Assumption about the minimal basis can be obtained with the 



help of Factorlnteger function from MATHEMATICA applied to the denominators. Such 
analysis allows to suggest the following basis: 



*W = (*■*• aFTT' 1 } < 310 > 



and it is easy to find the general expression for Eq. (13. 9ft : 



So, the contribution, proportional to C F nj is equal to 

7SU„; W = 16 Cmf^L - 1 Sl - | S 2 + I S 3 - i ^^ , (3.12) 

which is a full agreement with known result [281 ES]. Note, that constant is the same 
both for non-singlet anomalous dimension (13.71) and for non-singlet transversity anomalous 
dimension (J3.12J) . 

3.2 C F 2 nj contribution 

Now we proceed with the same method to reconstruct the rest parts. Again, we substitute 
N_ — > 1 and N + — > 1 into the result for the three-loop non-singlet anomalous dimension 
from Ref. [7] and obtain for C F 2 rir contribution 

-( 2 ) (m\ ian 2 f 23 55 c 5 Q L C L 2 c 20 c 4 c 
<W*0 = WCf nf KT6 - 24 Sl - 12 S * + S3+ 3 S * ~ Y S - 3 + 3 S - 4 

+ 4 9 4 9 + 4 °9 + 2 °9 8 9 4 9 

+ y S 2tl - ~S 2 , 2 - ^3,1 - l^a,! + (2 S t - ^ Cs) • (3-13) 



Subtracting the last equation at AT = 1, . . . , 15 from calculated values Eqs. ( 1A.lj) -( fA~T5|) 
we have found the following result: 

r , 2 f 56 4 8 5 206 14 313 761 3287 671 31789 



F f I 9 ' 3 ' 9 ' 9 ' 675 ' 45 ' 2205 ' 3780 ' 42525 ' 4725 ' 686070 ' 



6617 364087 



1171733 160379 1 , 
. 1, 10 n^ c . ^TT^TT^ ■ (3-14) 



62370 ' 12297285 ' 14189175 ' 8108100 J ' 
Again, analyzing denominators of obtained numbers we can write the following basis: 

« S- -J— — l ^ Sl 1 m^ 

c ^ n / \JV" JV + l' iV 2 ' (AT + 1) 2 ' AT' iV+lj ' { } 



With this basis it is easy to find the general expression for Eq. (I3.14p for even values of N 

^ 2 R )+ 2 (JV) = —C F \ , Sl - s (3.16) 

' TR > ns ' c F 2 V ; 3 F f N(N + l) v ; 

and for odd values of N 

~(2)- ,, n _ 16^2 ( Si 7 1 2 \ 

7tr^s,C7,V J~ 3 G ^ n /^iV(ivTI) 3iV(JV + l) (iV + l) 2 J- ld ° 

Note, that for (3 part our reduction rule gives correct result immediately 

^Z^JN) = C F \ (32ft - 24) Ca • (3.18) 

The final expression for C F 2 Tij contribution can be found from Eqs. (J3.13J) . (I3.16P and 

dSHD: 

3.3 C A C F rif contribution 

For C A C F rir contribution our reduction rule N_ — > 1 and N + — > 1 for the three-loop 
non-singlet anomalous dimension from Ref. [7j gives 

-( 2 ) (at\ inn r ( 5 209 q 167 q 7 c L 2 c 10 c 

lL!c A c F n f (N) = 16 C A C F n f { - - - —ft + —ft - -ft + -ft + - SL 8 

2 2 2 20 2 

~~ 0^-4 + ^1,3 — q^- 3 ' 1 + o'S'1,-3 — ^rft,-2 + q'S'2,-2 

-^3,1 + ^1,-2,1 - (2 ft - ~) Ca) (3.20) 

and subtracting the last equation at N = 1, . . . , 15 from calculated values Eqs. ( 1A.lj) - (1A.15j) 
we obtain 

J 22 35 31 7 371 293 2917 3877 47 1009 
A pUf I T ' 27 ' 54 ' 20 ' 1350 ' 1890 ' 17640 ' 45360 ' 420 ' 18900 ' 

74471 117259 340051 1470907 11840 l 
914760 ' 3243240 ' 5465460 ' 56756700 ' 2402400 J ' ^ ' 



With the same basis as in previous case ( 13. 15|> it is easy to find the following general 
expression for Eq. ( 13.2 lj) for even values of N 

lHt,c A c F n f (N) = ~-C A C F n f y N ^ +1) ~ YN(N + 1)) (3 ' 22) 
and for odd values of N 

7TR } nsr r n (N) =AC A C F n f ( — ^ + -—^ r - -7 ^ — ^ ) . (3.23) 

lTR,ns,C A C F n f K J A F * \N(N + I) 9 N(N + I) 3 (N + l) 2 J K ' 

The final expression for C A C F rir contribution can be found from Eqs. (I3.20p . ( I3.22p and 
fl3T23|) : 

4R%,C A C F n f (N) = lSc A C F n f (N) + l^c A C F n f (N) . (3.24) 



3.4 C F contribution 

For C F 3 contribution we again start with substitution N_ — > 1 and N + — > 1 into the result 
for the three-loop non-singlet anomalous dimension from Ref. [7] and obtain 

%!s 2) cv (N) = 16 C/ ( - j| + jjs 2 + ^3 - 4S 5 + 3S_ 2 + 95_ 4 + 12S_ 5 

-245_4,i - 4S_ 3 ,-2 + 65_ 3 ,i - 4S_ 3 ,2 - 125_ 2 ,_ 3 - 445 1> _ 4 - 6 Va 

— 1201,3 + 861,4 — 52o2,-3 — 6o2,-2 ~~ 662,2 — 4o2,3 — 12o3 j __2 ~~ 863, 2 

-4S 4i i + 24S_ 2 ,i,-2 + 64£ 1] _3,i - 165 li _ 2j _ 2 - 12S^ 2il + 8Si,_ 2 , 2 
+805 , 1 ,i,_3 + 165i, 2 ,_ 2 + 8Si j2j2 + 165i, 3 ,i + 56S 2 ,_ 2 ,i + 165 2)1] _ 2 

+852,1,2 + 85-2,2,1 + 85 3 ,i,i - 96S'i, 1 ,_2 1 i " (y + 12 ^- 2 ) Ca) • (3-25) 

Subtracting this equation at N = 1, . . . , 15 from calculated values Eqs. (1A.1[) - (1A.15[) we 
obtain 

„ 3 f244 > „ .10 A 1909 16 > 137437 48 _ 

CV i 136 Cs , 76 - 24 C 3 , 24 U , Cs , Ca , 

F I 3 S ' S ' 3 S ' 72 3 s ' 27000 5 S ' 

31939 12 113172503 36 1330189999 48 843934429 16 
2250 ~~ T ^ 3 ' ~ 27783000 ~ T ^ 3 ' 148176000 " 35 ^ 3 ' ~ 266716800 ~ y ^ 3 ' 

1950321041 8 2081091038207 24 2456339454271 48 
312558750 ~ 9 ^ 3 ' ~ 832031392500 ~ IT ^ 3 ' 532500091200 ~ 77 ^ 3 ' 
532264770130798741 144 729203803284910783 6 



263228107582440000 91 204732972564120000 13 

34163492121141653 6 > 1 , 

-7C3 • (3.26) 



20473297256412000 5 



Let's first try to reconstruct C3 part. It is clear from the last equation, that the possible 
basis has transcendentality level 1. However the basis from C F n? part f)3.10p does not 
work. Then, we slightly extended this basis 

Bc f s ^ = r 1 ' iv^r ' n ' w+i ' y+2/ (3 - 27) 

and we have found for even values of iV 

<2>+ ,„>_,„,- .,/ 1 1 \- %C F 3 ^ 



«Wc 8 W = 32 <>,'<, {—,-—2) = (jy-Dy+2) • < 3 ' 28 > 

while for odd values of N 

^;,, C ,.c,(^) = 288 C F % (I - ^ = |^| . (3,9) 

8 



For construction of the basis for the rational part we analyze again the denominators 
with Factorlnteger, for example: 

204732972564120000 = 2 6 x 3 6 x 5 4 x 7 4 x ll 3 x 13 3 . (3.30) 

So, we assume, that the basis for the reduced C F 3 part (I3.26P should contain the basis 
for C3CV 3 P ar t (13.271) . multiplied by harmonic sums S^N) up to weight 3. Moreover, 
we assume, that these harmonic sums should satisfy a generalized Gribov-Lipatov reci- 
procity [63 [68], which means, that up to weight 3 only the following combinations of the 
harmonic sums, called reciprocity- respecting harmonic sums [68, 69J g, should be consid- 
ered: 

I Si, S_2, S l5 S3, S 1 , SiS_2, 2S_2,i~ S_3j. (3.31) 

Thus, we stop with the following basis containing 18 terms: 



Bn 3 



/. 



11 11 



1 N{N + 1) ' (JV - l) 3 (iV + 2) 3 ' (JV - 1) 2 (JV + 2) 2 ' (JV - 1)(JV + 2) 

O-2 S_2 Si Si S_2Sl 



N(N + 1) ' (JV - 1)(JV + 2) ' N(N + 1) ' (JV - 1)(JV + 2) ' N{N + 1) 

C C C2 C2 Q3 C3 

D_2>Ji Oi £>i £>i D-l 



(JV - 1)(JV + 2) ' A^(A^ + 1) ' (JV - 1)(JV + 2) ' N(N + 1) ' (N - l)(N + 2) ' 

S3 S3 2S_2,1 — S_3 2S_2,1 — S_3 



iV(iV+l)'(iV-l)(iV + 2)' N(N + 1) , (N-l)(N + 2) 



(3.32) 



However, we have only 7 even values of anomalous dimension (I3.26p . from which we exclude 
N = 2 moment, because multigluon operator starts from N = 4. To find 18 coefficients 
from 6 known values we have used the same method, which we applied for the recon- 
struction of a general form of six-loop anomalous dimension of twist-3 operators in the 
maximally extended J\f = 4 SYM theory [39]. Major observation that the coefficients in 
any anomalous dimension are usually the rather simple numbers, that is the equation for 
the coefficients is a Diophantine equation. Moreover, a lot of harmonic sums, which can be 
written down for the possible basis, absent in final expression, i.e. its coefficients are zeros. 
There is a nice algorithm, which can help to solve such problem. LLL-algorithm [70J is re- 
alized in many computer algebra systems and its usage for our purpose can be found with 
the help of MATHEMATICA, where such algorithm is realized with function LatticeReduce. 
Firstly, we calculate the values of all 18 terms in the basis (I3.32p up to JV = 14. So, we 
have 6 equations in the linear system for 18 variables. We eliminate 5 variables and we 
remain with one equation on 13 variables. According to the realization of LLL-algorithno 
we add to the 14 x 14 unity matrix the column with numbers from the remaining equation 



Note, that the reciprocity-respecting harmonic sums are equivalent to the binomial sums [56U38J . 
2 See Application at http: //reference .wolfram. com/mathematica/ref /LatticeReduce .html 



of our system (13 + 1 numbers) and MATHEMATICA gives: 

fe.w = i«yk^+ „, fw', r ..» + s - 2 - es - 



\N(N+1) (N-l)(N + 2) N(N + 1) (N-l)(N + 2) 
GS-iS, | QS^ 2 S 1 6(2g_ 2;1 -g_ 3 ) \ (333) 



N(N + 1) (N-l)(N + 2) (N-l)(N + 2) 

Note, that LatticeReduce gives 14 possible solutions but all other solutions contain large 
integers and only few zeros. 

For odd values we should add to the basis few terms, which can be generated by 
Gribov-Lipatov reciprocity through a general equation [5"7] 155]: 

7 (N)=v(N + a"f(N)), (3.34) 

where a is some numerical factor and V(N) translated into the x-space is the reciprocity 

respecting function J67J EH]. This equation means, that the part of anomalous dimension in 

n-loop order is generated by the low order results. A new term proportional to 1/N/(N+ 1) 

appears at the first time namely for odd values of anomalous dimension at two-loop order 

(see the last line of Eq. (13. 6p ). This term will generate at third order the contribution, 

which looks like 

d ( S^N) \ 

So, we should add to the basis f !3.32|) the following terms with harmonic sums: 

f i>2 Oi Oi 



liV(iV + l) 'iV 2 '(iV + 1) 2 



(3.36) 



Indeed, for this extended basis LLL-algorithm gives the following general expression for odd 
numbers from Eq. (13. 26ft : 



TR,ns,cvv / t ViV(Ar + i) (N-l)(N + 2) {N + l) 2 N 2 {N + 1) 2 

45i 55i 2^1 2S 2 S_ 2 

+ T77TT 7T ~ TTT 7T7TT ^ + T77TT 7T + 



(iV + 1) 2 N(N + 1) (JV-l)(iV + 2) iV(iV + l) (JV-l)(JV + 2) 

2o_2i- , i loO_2 1 - , l 18(2o„2 1 — <J— 3 



(N-l)(N + 2) N(N + 1) iV(iV + l) 



(3.37) 



The final expression for C F 3 contribution can be found from Eqs. f !3.25j> . f !3.28j> . (I3.29J) . 

( 13T33J) and (EOTj) as 



7?itc a (N) = 7^3 (N) + 7^1.^3 (AT) + tSI ,, .,, (N) ■ (3.38) 
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3.5 C A C F contribution 

Substitute N_ — > 1 and N + — > 1 into the result for non-singlet anomalous dimension from 
Ref. [7] we obtain for C A C F 2 contribution 



'ns,c A c F 2 v ; a f y 64 24 2 6 2 



7 n 7r r »W = 16C A Cy - — + -rrS a - —S 3 + -5 4 + 105 5 - -5_ 



2 



1 34 8Q 31 

+ ~^~ "-3 (T^~ 4 ~ 1^5-5 + 20S'_4 ) i — 2S'_3 ) _2 —^-3,1 + 2o_3 : 2 

iin , fiQ , ^ 31 268 ' 134 
+106_2,_3 - 66_ 2 ,-2 + 345i _ 4 + — 6i _ 3 —S 1 _ 2 g-^1,2 

+ — 5*1,3 — 22S'i i 4 + 425*2, -3 + "TT 52,-2 — ^2,1 + ~ 52,2 ~ 1052,3 

6 6 » 6 

4 
+105 3i _ 2 + tt53,i — 25 3)2 — 2S , 4> i + 85_2,-2,i — 285_ 2 ,i,_2 — 645i,_ 3) i 

+32S'i i _2,-2 + — 5i, -2,1 — 451,-2,2 — 7201,1,-3 + 165i,i,3 ~~ 801,2,-2 

-165i, 3i i - 605 2 ,- 2 ,i - 85 2il ,- 2 + H25 1 , 1 ,_2,i + \ (y + 125_ 2 ) Cs) • (3.39) 



Subtracting this equation at N = 1, . . . , 15 from calculated values Eqs. (lA.lj) - (1A.15j) we 
obtain 



C A C r > { - ifi + 220C3 , -f + 36 C 3 , 1 + 36 C 3 , -*£ + 8< 3 



, j 1424 : oonA 283 , Q ^ A 5 QfiA 4477 , cA 94739 72( 3 

" '' ' 18000 + ~1T 



36817 , I8C3 28787477 , 54C 3 3040082519 , 72( 3 2007772853 24( 3 



2250 5 ' 6174000 7 ' 296352000 35 ' 533433600 
11821223651 4( 3 20222794198253 , 36Cs 27807644469347 , 72( 3 



1666980000 3 ' 6656251140000 11 ' 5325000912000 77 ' 

1310584038208950541 216C 3 1646199071310197233 9( 3 
526456215164880000 + 91 ' ~ 409465945128240000 + "l3 ' 

8492728741110641 % 3 ] , 

+ — }. 3.40 

4094659451282400 5 J v ; 

Using the same basis (I3.32p as in previous subsection we have found with LLL-algorithm 
the following general expression for even numbers from Eq. f )3.40p 



^Zs,C A cA^ = 8C ^ 2 { 



1035! I8S1 35_ 5 



6N(N+1) (N-l)(N + 2) N(N + 1) 



185_ 2 185_ 2 5i : 65_ 2 5i 18(25_ 2 ,i — 5_ 3 ) , , , 



(N-l)(N + 2) (N-l)(N + 2) N(N + 1) (N-l)(N + 2) 
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For odd numbers from Eq. (13.401) we extend basis (13.321) with few additional terms ( 13.36 
and obtain 



^tL,c a c^ N ) = 8C aC f 2 ( 



31 209 3 451 

+ 



3(iV+l) 2 9N(N + 1) (JV-l)(JV + 2) (N + l) 



2 



2b\ 5oi>x Zbi 2o2 



N 2 (N + 1) 2 QN(N + 1) (N-l)(N + 2) N(N+1) 

35„ 2 + 25Ugi _ 54S_ 2 S 1 + 54(2S_ 2|1 - S_ 3 ) , ( , 



(N-l)(N + 2) (JV-l)(JV + 2) N(N + 1) N(N + 1) 

For ( 3 part the results are the same as in C F 3 case in Eqs. (I3.28P and (I3.29|) up to 
common factor (—3/2). The final expression for C A C F 2 contribution can be found from 
Eqs. ( 13391 . (JS3I]) , ( M) , dOSD and ( 137291 as 

^ 2 (iV)=7 (2) 2 (iV)+7 (2)± a (jyn--S^7 (2)± 3 (N) (3 43) 

'TR,ns,C A C F 2 v J V 'ns,C j4 C F 2 v 1 V ^ 'TR,ns,C* A C* F n ^ 2 C 3 r TR,ns,C F 3 C3 ^ ' ' \°-^°) 

3.6 C A 2 C F contribution 

For the last C A 2 C F contribution our reduction rule N_ — > 1 and N + — > 1 for the result of 
the three-loop non-singlet anomalous dimension from Ref. [7] gives 

^ 2 )+ (N\ i fi r V f 1657 u 245 ? 1043 389 31 

7 ns ,c/c F (^) = 16 C A C^— + —ft - —S 2 + —S 3 - -S 4 - 2S 5 

+ls_ 2 - ^SL 3 + ^S- 4 + 2S_ 5 - 4S_ 4 ,i + 25_ 3 ,-2 + ^S-3,1 
2 9 6 3 

— 2S'_2,-3 + 3S'_2,-2 — 601,-4 — "1,-3 H ^-"1,-2 ^"*S'l,3 + 65*1,4 

3 y 2 

— 85*2 3 — 1S2 _2 + 4S"2 3 — 2S*3 2 + "7" "3,1 — 2S*4 i — 45_2 _2,1 

3 6 

22 
+85'_2,i _2 + 1651,-3,1 _ 125'i,_2,-2 — ttSi,-2,i + 165i,i,-3 — 851,1,3 

+85-1,3,1 + 1632,-2,! - 325i,i,_ 2 ,i ~\(y + US - 2 ) Cs ) • (3 ' 44) 

Subtracting this equation at N = 1, . . . , 15 from calculated values Eqs. (1A.1[) - (1A.15[) we 
obtain 

A F l 3 27 3 120 3 5400 5 

134749 _ 6Cs _ 846463 _ I8C3 14190767 24( 3 10288711 __ 8(3 
37800 ~ IT '~ 396900 T~ ' 6350400 ~ ~35~ ' ~ 6350400 ~~ ~5~ ' 
1803287 4C 3 972218609 12C 3 85376995477 24( 3 
1166400 ~~ ~Y ' ~ 768398400 ~ ~TT ' 74918844000 ~~ ~W ' 

593181857171 72( 3 4789546276259 3C 3 7574874351263 3C 3 
"584366983200 ~ IJT ' 5454091843200 ~~ "l3 ' ~ 9090153072000 ~ ~T 
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(3.45) 



Using the same basis (I3.32p as in subsection 13.41 we have found with LLL-algorithm the 
following general expression for even numbers from Eq. ( 13.451) 

^ + (N) = C 2 C ( 121 I 56Sl I MSl 

7 T R,ns,cVC F W - W ^F y giv(iv + x) I" 3iy(iv + 1) + (JV - 1) (JV + 2) 

4g- 2 24S_ 2 245^1 _ 24(2S_ 2 , 1 -g_ 3 ) \ ( , 

N(N + 1) (N - 1)(N + 2) (N - 1)(N + 2) (N - 1)(N + 2) / ' [ ' 

while for odd numbers from Eq. f )3.45[) we obtain with extended basis from Eqs. ( I3.32p 
and (ESHD 

^ 2) ~ (N)-C 2 C(- U 25 4 



3(iV + l) 2 N(N + 1) (N-l)(N + 2) 
4Si 4S_ 2 72g- 2 g! | 72(2g_ 2il - ff_ 3 ) , ( , 



iV(iV + l) (JV-l)(JV + 2) iV(iV+l) iV(iV + l) 

For £ 3 part the results are the same as for C F 3 contribution in Eqs. ( I3.28P and ( I3.29P 
up to common factor 1/2. The final expression for C A 2 C F contribution can be found from 
Eqs. (EBD, (J33SD, (E3ZD, (J52BD and (ESSD as 

4 Results in Mellin space 

Here we present the final results for the three-loop anomalous dimensions of the flavour 
non-singlet transversity operator (12. ip . For representation of our results we use the same 
notation as in Ref. [7] 

N±^ = ^(iV±l), N ±i ^ = ^(iV±z) (4.1) 

for arguments shifted by ±1 or a larger integer i. In this notation the one-loop (LO) 
anomalous dimension [191 EQl EH [181 [22] reads 



1 { tLs(N) = C F (AS 1 - 3) (4.2) 

and the corresponding second-order (NLO) non-singlet quantities [2H ESJ [26] are given by 

1 { tIU N ) = 4 C aC f (-^ + Y Sl ~Y S2 + 2Ss ~ 2S - 3 + 45l <- 2 
+4C F r V (^-y5 1 + ^ 2 )+4CV 2 (35 2 -| + 45_ 3 -85 1 ,_ 2 - \S L2 IS 2A ). ,1.3) 

,(1)- fAT\— n ,W+ fAT\ I Q^ (n C A 



1$£JN) = 1tAU N ) +8C F (C F -^j (N_ + N + - 2)5 t . (4.4) 
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The three-loop (NNLO, N 2 LO) contribution to the anomalous dimension 7TRns(-^0 f° r 
even moments reads from Eqs. (JS3D, (13~T3l . (13T2DJ) . (13T25D . (I3T39J) . (13^4j) as 



+16C„ 2 n 



r - ; U( N + - 1)^2 - J(N_ + N + - 2)5i,i 



+ 16C A C F n^ - -(N+ - 1)6 2 + (N_ + N + - 2) 

C, 



11 9 + N 

IT 1 + l^ 1 ' 1 



+32C F (C F - C A ) (C F - -^J (N_ 2 - N_ - N + + N +2 )6iC 3 



H-16C. 2 CU (N+ - T 



— i>2 — — 63 
3 2 A 



-(N_-N + )6 2 ,i + (N +2 -N 
+ (N_ 2 -N_-N+ + N +2 
+16C A C F 2 (3(N_ - N + )6 2 ,i + (N+ - 1 



+ (N_ + N + - 2) 

5*2,-2 — ^3,1 



— S 10 5 -^6 



"31,-2,1 — "1,1,-2 

139 



"2 + — 63 + 3o4 + 3o2 : -2 — 3o3 : i 



+3(N +2 - N 

— 3Si 2 ,i ~~ 3oi 1 _2 



"3,1 — 5*2,-2 



(N_ + N + - 2 



12 2 

139 



61 i — — 61 _o ~T 36 1 _3 

12 ' 2 ' ' 



3(N_ 2 -N_-N + + N +2 ; 



+16C F 3 ( - 2(N_ - N + )6 2 ,i + 2(N +2 - N_ 



6- 



2,-2 



5i,i, 

5*3,1 



5Y 



-2,1 



+ (N+ - 1 
+ (N_ + N + -2) 
+2(N_ 2 - N_ - N + + N +2 
For odd moments we have 



62 — "3 — 664 — 662,-2 + 663,1 

8 "5*1,1 + 5i 5 _2 — 65i 5 _3 + 66i 5 _2,i + 601^,-2 

"31,-2,1 ~~ "31,1,-2 



(4.5) 



4tns(N) = 4tU N ) + WC F n f (C F - -^)(N_ + N + - 2) 



C. 



7 2 

~ 9 Sl + 3 s1 ' 1 



+ 



C_ 



+32C F (C F - C A )(C F - -^)(9(N_ + N + - 2)6 X - (N_ 2 - N_ - N + + N +2 )6i)C 3 



+16C A C F (c F - ^f)((N + - 1)(26 2 - 5 3 - 1862,-2 + 185 3 ,i 



+ (N_ - N_ 



762 — 262,1 



+ (N_ 



N_ 



-63 + 262,-2 — 263,1 



+ (N_ + N+-2; 



179 17 

——61 + 61,-2 — 5~5i,i — 186i,_2,i + 186i,i, 

lo o 



14 



0.25 
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Figure 1: The perturbative expansion of the anomalous dimension 7TRns(-^0 f° r f° ur 
flavours at a s = 0.2. In the right part the difference between our three-loop result and the 
corresponding result from Ref. [7J is shown. 



+ (N_ 2 - N_ - N+ + N 
+16C/(C F -^)((N + - 



+2) 



-S\-2 + 251,-2,1 — 25i ) i : _2 

19 
- 45 2 + — S 3 + 6S 4 + 245 2 ,_ 2 - 45 2l i - 24^3,1 
o 



+-(N_-N + )S 2 , 1 + (N +2 -N_ 



+ (N_ 2 - N_ 
+ (N_ + N + -2 
— 2oi 9 + 12iSi _2,i 



N + + N +2 
13 9 
245 1 ,i_ 2 



12 8 8 

Q Q q 1 C 

o^3 — o»J4 — o^2,-2 + o«->3,l 

o^l,-3 + o'S'l-2 — o'S'l-2,1 + TT'S'l,! -2 



0^2 — -^Ss + 661,-3 



5 9 + 23 9 



(4.6) 



The results ( 14. 2p . H4.3|) and (14. 5 p for 7tr, ns (iV) are collected in Fig. [TJfor four active 
flavours and a typical value a s = 0.2 for the strong coupling constant as in Ref. [7J. 
Numerically, the colour factors take the values C F = 4/3, C A = 3. Moreover, in the right 
part we show the difference between 7TR ns (iV) and y£(N) from Ref. [7J and one can see, 
that the difference between NLO and NNLO results is very small. 
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5 Results in x-space 

The non-singlet anomalous dimensions 7TRns(^0 are related by the Mellin transformation 



7 ( t£ s (A0 



dxx 



N-l p(n)± 



P TRnsW 



(5.1) 



to the splitting functions -P^Rns^)- So, the splitting functions P^R ns {x) can be obtained 
from the iV-space results of the previous section by an inverse Mellin transformation, 
which expresses these functions in terms of harmonic polylogarithms [TTJ [721 E3]- The 
inverse Mellin transformation can be performed by a completely algebraic procedure with 
HARMPOL package [5S1 [7J for FORM [SO]. 

Our notation for the harmonic polylogarithms H mi) .„ >mw (x), rrij = 0, ±1 follows Ref. [73] 
to which the reader is referred for a detailed discussion. The lowest-weight (w = 1) 
functions H m (x) are given by 



Ho(x)=\iax, H±i(x) = =F hi(l =F x) . 
The higher-weight (w > 2) functions are recursively defined as 



(5.2) 



— In™ x . 



£lm\,...,m w \%) 



W ! 



px 

j dzf mi (z)H, : 
Jo 



m2,...,m w 



(*) 



\im x ,...,m w = 0, 
else 



with 



/o(x) 



f±i{x) 



x 1 =f x 

A useful short-hand notation is 

Ho, . . . , o, ±1, 0, . . . , 0, ±1, ... (%) — #±(m+l), ±(n+l), 
m n 

The one-loop non-singlet splitting function [TjJJ [2H1 1211 dSl [22] can be written as 



(x) 



°TR,ns( a 



where 



C F (2p qq (x) + 35(l-x)), 
2x 



Pm( 



x 



1 — x 



(5.3) 

(5.4) 
(5.5) 

(5.6) 
(5.7) 



Here and below we suppress the argument x of the polylogarithms and all divergences for 
x — > 1 are understood in the sense of +- distributions. 

The two-loop non-singlet splitting functions [2U [251 I2S] are given by 



P TR.ns( X ) 



+4C A C F [p qq (x) 



AC F n f p qq (x) 
67 



IS 



-C 2 



5 

9 3 
11 



iH (1 + *(i-,)(-L + !fc) 



6 



H + H ,o 



+ Pqq(- 



-.c 



C2 + 2H_i >0 — H ,c 



16 



+AC F 2 [2p m (x 
+45(1 -x) 



Hi,o — tHq + H 2 



SPqql-^) 



^,^V(g + fc 2 -3<3)+CV 2/: ' 



(2 + 2H_ li0 — H ,o 
- 3C 2 + 6C3 



(5.* 



°TR,ns( 2 ') 



+16C F (C F -^)( Pqq (-x) 



-* TR,ns( X J 



C2 + 2H_i )0 — H ,i 



1, 



(5.9) 



The three-loop splitting function for the evolution corresponding to the anomalous 

/Q\_|_ 

dimension 7 T R ns (iV) from Eq. (J4.5P reads 

1 

27 



-* TR,ns( X ) 



16C>/(— (l-x) + — + /> m {.r) 



54 a; 54 ° 18 ' 



.r 



+16C A C F n f ^-—n + (l-x 
+ 7C2H0 - -H 3 



18 12 



Pqq(-Z) 



1 5 2 

9^ 3 ~~ cp 2 — 0C2H-1 



6 

+Pqq(^) 



10 H + 2 H + 5 H X H + X H 

~5~ n -i,o + o n -i,2 t q n o,o — o • n —i,o,o ~r -no,o,o 



+ lQC F 2 n 



0—0 o-^-2,0 

Cs " 2 Cs - 216 + (s C2 " Tos) H ° " 6 H3 " 6 H °'° 

C2 — C3 + 77C2H-1 — 77C2H0 + - H 3 



9 



ttHo,o,o + 7^1,0,0 



X , 



-H r 



F /v 3 ~» 3 



2 20 4 

+ 0^-2,0 + -q-H_ 1)0 - -H_i )2 



1 -x)Hi +p m (-x 
10 



9 



^-H 0> o + 7tH_i )0 ,0 



-H, 



H-PqqO) 



5 55 /l 

3^ 3 48 + V3^ 2 



I Hn — H? — — Hq 

24/ 9 3 



0,0,0 



X H 10 H 2 H 



X H 4 H 

— o 110,0,0 — o-H-1,0,0 



16C/ H ,o - 3C 3 - 8H + 6C2H0 - 6H 3 - 6H_ 2i0 + 6H , ,c 



+x* 



C3 — 2C2H0 — 2C2H1 + 2H 3 + 2H_ 2 ,o 



+ (l-x) 



6C3 - 2C2 - 9C 2 H - 8Hi 



-3C2H1 + 6H 3 — H oc 



+ [ - + r 

X 



C2H_i + C2H1 — 2H_ 12 — 2H_i 1 



9Cs - 2C2 - 9C 2 H_! + 8H + 3C 2 H + 2H 2 + 6H_ 2>0 + H_ 1>0 + 6H_ li2 



+ (l + x 

— 6H_i _i )0 + 6H_i )0 ,o — 6H ,o,o 



Pqq(-X) 



-A2 



9 



-Cs + 32C 2 H_ 2 + 6C 2 H 2 + 3H 3 



+ (6(2 + 36C3)H_! + (-13C 3 - \ - 7^2) H + 12H 4 - 6H_ 3 ,o + 3H_ 2 , - 28H_ 2 , 2 

-48C 2 H_ 1 ,_ 1 + 40C 2 H_ 1>0 - 6H_ li2 - 32H_ li3 - 14C 2 H , + 2H 3 , + 8H_ 2 ,_ 1)0 

9 
— 26H_ 2j0 ,o + 8H_i 7 _2,o + 48H_i ) _i )2 + 3H_i i0 ,o — 4H_i, 2)0 — -H 0; o,o + 40H_i i _i >0 ,o 

-22H_ li0 , ,o + 6H 0i0 ,o,o] + P qq (x) [^2 2 + 6C 2 H_ 2 + (Cs ~ |Ca - ^) H + ^3^ 

13 
+2H4 — 2H_3 5 o + —Hq^ + 8Hi,3 ~~ 3H 2j o + 4H 2 , 2 + 4H3,o + 4143,1 + 12H_2,-i,o 



17 



— 6H_ 2 ,o,o + 8Hi 5 _ 2 ,o — 6Hi ;0 ,o + 4Hi ;2 ,o + 2H 2 ,o,o + 4H 2 ,i,o — 2H 0) o,o,o — 4Hi )0 ,o,o 

+16C A 2 C F { - ^H + 1h , + x 2 [ - |Cs - C2H0 - C 2 Hx + H 3 + H_ 2 ,o^ 

., 1 , sfl21 1 u 10„ 1 



■ x) 



+ f - + x 2 
\x 



- TC2 + y H + H 2 + - 

n 



nH-1,0 



(l-x) 



~C 2 H_i + 2C2H1 - H_i ;2 - H_i _i ;0 + p qq ( 



121 V 10 H X H 
"72" " 4C2 - yHi - -H 0i0 

[67 . 2 11 
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o 
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9 
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31 



— ^-H_i j2 - 8H_ 13 
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+ Pq.q( X ) 

Ho 

—I 

12 



11 31 
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2 s 5 18 48 12 
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/ A 31 1043\ 
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/- /- /- TT /"- TT 

4C2 - -C3 + T^H- 1 ~ 2^ 2 ° 
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Figure 2: The three-loop splitting function ixRnsO*') f° r different numbers of active quarks, 
multiplied by (1 — x) as in Ref. [7J. In the right part the difference between our three-loop 
result and the corresponding result from Ref. [7j is shown. 
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The function P^ ns (x) from Eq. (15.101) is shown in Fig. [2] together with the difference 
between our result and the result for P^ + (x) from Ref. [7J. The numerical evaluations of 

/o\_i_ 

the harmonic polylogarithm, entering into expression for -RpRns^) fr° m Eq. ( I5.10p . were 
performed with HPL package [75j [76] for MATHEMATICA. 

The three-loop splitting function for the evolution corresponding to the anomalous 
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dimension 7TRns(-^0 fr° m Eq. (14. 6 p is given by 
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6 Conclusion 

We have found the complete third-order contributions to the splitting functions governing 
the evolution of flavour non-singlet transversity parton distribution in perturbative QCD. 
Firstly, we calculate the fifteen moments of the anomalous dimension of the non-singlet 
transversity operator with the help of MINCER package [42] for FORM [5D]. As in our 
previous paper [3H] for the reconstruction of the complete anomalous dimension for arbi- 
trary moment N from the know fixed values we have used LLL-algorithm [70], that allow 
to find coefficients in the ansatz from the harmonic sums, which form the suitable ba- 
sis. We have assumed, that the most complicated part of anomalous dimension could be 
obtained by replacement N_ —$■ 1 and N + — > 1 in the known result for the non-singlet 
anomalous dimension from Ref . [7J . We used the same procedure in Ref . [36J together with 
the maximal transcendentality principle for the computation of the three-loop universal 
anomalous dimension in the maximally extended M = 4 SYM theory In present paper, 
we have subtracted this part from the known fifteen moments of anomalous dimension, 
listed in Appendix. The obtained numbers look much more simple with compare to the 
original numbers from Eqs. (1A.1[) - (1A.15[) and we can expect that the basis will reduced 
considerably. Apply LLL-algorithm to the minimal basis and the highest even or odd ob- 
tained numbers we have reconstructed the complete three-loop anomalous dimension of the 
flavour non-singlet transversity operator for arbitrary N. From obtained expressions the 
x-space splitting functions can be obtained by a standard Mellin inversion with HARM- 
POL package [23 [75] for the symbolic manipulation program FORM [5U] and we present 
the corresponding results in section [5J 

For possible test of our result we give below an explicit expression for the three-loop 
anomalous dimension 7TR ns (^0 from Eq. (14.31) at iV = 16: 



(2) . , 9588431966492629867^ 2 

'TR, ns v i 3579902263121184000 F f 



( 26209369318544258455663 3792038 > 
-C A C F n f —————— — — — + — — Cs 



f V 1002372633673931520000 45045 
/214170060836708111694659 379203^ 
%f V 2810574247360239360000 45045 



5187682484935485448445159707 100658462371 
303845560733620756730880000 ~ 2705402700 ^ 3 
2 /24000075522445759616746817 100658462371 
V 590220591945650265600000 + 1803601800 ^' 



, ,' 1448460151630990534353391 100658462371 



+ A F \ 5613286748574016512000 + 5410805400 ^) ' ^ ' ' 

The result of the recent calculations for iV = 16 Mellin moment of the three-loop 
anomalous dimension of the non-singlet transversity operator from Ref. J4T] coincides with 
the above prediction and can serve as a direct confirmation of the general result presented 
in this paper, at least for even values of N . 
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FORM file of our results can be obtained from the preprint server http : //arXiv . org 



by downloading the source. Furthermore it is available from author upon request. 
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Appendix: The 3— loop Anomalous Dimensions 

7SU(1) = ~\ C F n?T* - n f T F (c a C f (^ + I6C3) + C/ (^ + 16& 

+c/ (f _ 64C3 ) - c A cS («f» - n%) + c/c F (i*» - 4o C 3 ) , (a,; 

tSU(2) = -p C F n f 2 T F 2 - n f T F (c F 2 (30 - 48C 3 ) - C A C F (^ + 48, 



+ ^" ( - / f -t^a^f "^ Tns - ^ F ' ( A -2) 

^(3) - -^ w - ^ K (H + *?) + <v (i||f + *?)) 

"^ ( W - C,f ) - c A c, (^1 + 32C3 ) + csc f (^ + 6 f) . ,a,) 

(2) , A . 7361^ 2 ^ 2 ^ /^ ^ /199723 A 256\ ,_, 2 /66443 A 256\ 
7 W (4) = CpntTl - n f T F C A C F + Cs + C F 2 [ + (3 

/TR,nsV ) 1215 F f F f F^ A F^ ^gg S3 g J F y Q72 S3 g y 

-,/ (^ - 2SC) - ^ (i|i|i + <*) + ,/C (= + .0) , (A.4) 

/ 49495163 . 1472\ „ , / 15920231 . 2356 



-%2><V — — — - c 3 — - C 



/i/V ; V 607500 S3 15 J ^ F V 759375 S3 75 ) 
2 /65092847 1178\ „ 2 ~ /79012363 u 1178 



~ CaCV I 2430000 + C3 ^5~J + ^ ° F { 486000 + Cs ^5~ } (A/j) 

(2) , R , 1604879 22 /1841332 544 

7tr „, 6 = Crp-nfTp — n f TpC A C P h C3 

mi,nsi ) 212625 F f F f F A F \ 42525 s 5 

2 f 18622301 544 \ 3 / 4075862 832 

~ Uf F F V 202500 ~ ^ 3 ~) ~ F \ 253125 ~ ^ 3 l25~ 

A F V 67500 S 25 / A F V 63000 S 25 ' 
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(2) 84250571 22 / 1844723441 4128 \ 

7TR ' ns(7) " "10418625 °^ ^ ~ n /^ C ^ (, 41674500 + G ^5"J 

2 / 49282560541 _ 4128\ _ 3 /232913441459 _ 42208\ 
" n/ F F V 486202500 " 135/ ~ F \ 11344725000 ~ ^ 3 1225 ) 
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~ Ca ° f { 648270000 + C3 l225~J + ^ Cf { 27783000 + C3 7225~J ' (A ' 7) 
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~ Hf F F V 7779240000 ~ ^ 3 35 ) ~ F \ 136136700000 ~ ^ 3 1225 ) 



—C A C F 
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+C A 2 C F [ h C3 , (A.10) 
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(2) = _™™ 
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— C C h c 
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^ „ „ ( 1063558335824080583 , 721 1704 \ 

-n f T F C A C F + Cs 

f F A F V 20938599466785000 s 45045 / 
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/ F F V S7«^n44.ci9^nnn74snnnn ^ 3 



V 87830445230007480000 ~ 45045 
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